Exploring Zingiber officinale Essential Oils Compounds as Selective Inhibitors of Orexin Receptors: A Comprehensive Molecular Docking and Dynamics Study
Raharjo SJ1*, Sari DRT2, Wijayanti ED3, Maryanty Y4
1Department of Pharmacy Analyst and Food, Health Polytechnique of Putra Indonesia Malang, Indonesia.
2Department of Pharmacy, Faculty of Health Science, Ibrahimy University, Situbondo, Indonesia.
3Department of Pharmacy, Health of Polytechnique of Putra Indonesia Malang, Indonesia.
4Department of Chemical Engineering, State Polytechnic of Malang, Indonesia.
*Corresponding Author E-mail: sentotjoko@poltekkespim.ac.id
ABSTRACT:
Insomnia is the second most prevalent mental disorder and is associated with various health risks. Ginger (Zingiber officinale) essential oil is widely recognized for its diverse therapeutic properties. This study aimed to explore the potential of compounds derived from ginger essential oil as therapeutic agents for insomnia by targeting orexin-1 (OX1R) and orexin-2 (OX2R) receptors through molecular docking and molecular dynamics simulations. Forty compounds identified from ginger essential oil via GC-MS analysis were retrieved from the PubChem database, while the crystal structures of OX1R (PDB ID: 4zj8) and OX2R (PDB ID: 7xrr) were obtained from the Protein Data Bank. These compounds were docked to the active sites of both orexin receptors, and their binding interactions were compared with those of native ligands. Molecular dynamics simulations were performed using OpenMM on Google Colab to assess the binding stability of the compounds using ligands with the best binding affinity. Docking results showed that ginger essential oil compounds bind to the same sites as native ligands on both receptors, exhibiting relatively higher binding affinity and selectivity to OX2R, with alpha-zingiberene identified as the best ligand. MMGBSA and MMPBSA binding energy calculations from molecular dynamics simulations showed ΔG values of -8.9200 ± 2.3693 kcal/mol and 0.0226 ± 3.1674 kcal/mol for the alpha-zingiberene–OX2R complex, respectively. In conclusion, molecular docking showed that ginger essential oil compounds generally exhibit selectivity for OX2R. Molecular dynamics simulations also showed that alpha-zingiberene forms a stronger and more stable interaction with OX2R.
KEYWORDS: Insomnia, Essential oil, Ginger, Orexin receptor, Alpha-zingiberene.
INTRODUCTION:
The pace of life has accelerated considerably, affecting people’s lifestyles and daily habits. Nowadays, many individuals dedicate more time to work rather than rest1. This shift has become a major factor contributing to a restless lifestyle, which in turn leads to stress and disturbances in the circadian rhythm2. Excessive activation of the stress system results in a hyperarousal disorder known as insomnia. This condition is marked by difficulties falling asleep, impaired daytime functioning, fatigue, memory problems, and a persistent cycle of anxiety. Current insomnia treatments primarily focus on cognitive behavioral therapy and pharmacological approaches3. However, some studies have reported adverse effects associated with these treatments, including dependence, tolerance, and drowsiness4,5.
The sleep-wake cycle is primarily regulated by a pair of neuropeptides known as hypocretins, which are typically produced by orexin-containing neurons located in the lateral hypothalamus (LH). These neuropeptides interact with two types of receptors, orexin-1 (OX-1) and orexin-2 (OX-2), which are connected by small peptides and distributed throughout both the central and peripheral nervous systems6,7. When orexins bind to their specific G-protein-coupled receptors (GPCRs), they trigger downstream signaling pathways that carry out various physiological roles, including the promotion and maintenance of wakefulness8,9. Research conducted on freely moving cats revealed that OX-1 levels were higher during active waking periods in the daytime compared to slow-wave sleep (SWS), and these levels were also elevated during rapid eye movement (REM) sleep relative to SWS10. Another study found that OX-1 concentrations in cerebrospinal fluid (CSF) were low immediately after waking but gradually increased throughout the day’s activity and then declined slowly at night11. However, orexin expression is increased during stress, which can disrupt the circadian rhythm12. The Food and Drug Administration (FDA) has approved several dual orexin receptor antagonists (DORAs), including Suvorexant, Lemborexant, and Daridorexant, for the treatment of insomnia9. Despite the positive therapeutic effects of these drugs, frequent use has been associated with adverse side effects such as drowsiness, sleep paralysis, headaches, dry mouth, hallucinations, and complex sleep behaviors13. Therefore, exploring new potential orexin antagonists derived from bioactive compounds could help minimize these drug-related side effects.
Zingiber officinale, commonly known as ginger, is a perennial herbaceous plant renowned for its diverse medicinal properties. The bioactive compounds found in ginger exhibit strong anti-inflammatory and antibacterial effects14-16. More recently, ginger’s potential effects on the central nervous system have been revealed17,18. Its volatile compounds interact with amino butyrate aminotransferase, acting as substitute GABA molecules19,20. Specifically, 6-gingerol, a major constituent of ginger, has been shown to effectively reduce seizure-like paroxysmal behavior in zebrafish larvae. This compound is also believed to help restore the balance between GABA and glutamate (GLU) in the epileptic brain19. Additionally, ginger’s bioactive components are known to cross the blood-brain barrier18. The neurophysiological activities of ginger have prompted further investigation into its effects on various neural mechanisms. This study aims to evaluate the therapeutic potential of ginger essential oil compounds for treating insomnia by targeting orexin-1 (OX1R) and orexin-2 (OX2R) receptors using molecular docking and molecular dynamics simulations.
MATERIALS AND METHODS:
Preparation of Ginger Essential Oil:
The rhizome of ginger (Zingiber officinale) was first processed into a simple powder. Essential oil was then obtained from 250 grams of this powdered ginger through steam distillation. The oil was extracted using n-hexane and subsequently purified by evaporating the n-hexane solvent21.
Gas chromatography and mass spectrometry of Ginger Essential Oils:
Ginger essential oil was analyzed using gas chromatography-mass spectrometry (GC-MS QP2010 Plus) equipped with a non-polar/polar stationary phase column. The oil sample was injected into an Agilent 19091S-433 HP-5MS column. The column temperature was initially set and then increased gradually by 10°C every 28 minutes until it reached a final temperature of 288°C. The chemical components were identified through mass spectrometry (MS).
Receptors and Ligands Retrieval:
The interaction between Zingiber officinale volatile compounds and orexin receptors was investigated through molecular docking and molecular dynamics simulations. The structures of orexin-1 receptor (OX1R, PDB ID: 4zj8)22 and orexin-2 receptor (OX2R, PDB ID: 7xrr)23 were retrieved from the RCSB Protein Data Bank. Preparation of the receptors involved removing water molecules and native ligands, followed by the addition of hydrogen atoms using AutoDock Tools24,25. The 3D-structures of 30 selected ligands were downloaded from PubChem in SDF format, converted to PDB files, and subjected to energy minimization using Avogadro ver. 1.2.0.
Molecular Docking and Dynamic:
Molecular docking was conducted using AutoDock Vina version 1.2.3. The docking of ligands to the receptors was performed within a grid box optimized automatically by the software. This grid box was validated by re-docking the native ligand for each receptor to ensure accuracy. The resulting 2D and 3D docking poses were visualized using Discovery Studio26-29. To further investigate the binding interactions of the best-scoring ligand, alpha-zingiberene (CID92776), with orexin receptors OX1R (PDB ID: 4zj8) and OX2R (PDB ID: 7xrr). Molecular dynamics (MD) simulations were performed using OpenMM on the Google Colab platform. For ligand preparation, the 3D structure of alpha-zingiberene was obtained from the PubChem Open Chemistry Database. The SMILES representations of three different compounds were converted into PDB files using LigParGen, a web-based tool that generates force field parameters for organic molecules and ligands30–32. Before running the MD simulations, the necessary software libraries and packages were installed, including Anaconda, OpenMM, PyTraj, py3Dmol, ProLIF, NumPy, Matplotlib, and AmberTools. The ligand and protein PDB files were uploaded and managed on Google Drive and Google Colab. The MD simulations followed the protocol outlined by Pablo R. Arantes et al. Protein topology was generated using the FF19SB force field, with TIP3P water molecules. The system was neutralized with NaCl ions at a concentration of 0.15 M using AMBER tleap. Ligand topology was generated using the GAFF2 force field. The equilibration phase included 1,000 steps of energy minimization, followed by a 5-nanosecond simulation with a 2-femtosecond integration time step, at 310°K temperature and 1 bar pressure. Trajectory and log files were recorded every 10 picoseconds. Production MD simulations were then run for 10 nanoseconds with similar parameters: a 2-femtosecond time step, 310°K temperature, 1 bar pressure, and data recorded every 10 picoseconds. Binding free energy calculations were performed using both the MM-GBSA and MM-PBSA methods for comparison. The GB/SA parameters included the OBC model with igb=2 and a salt concentration of 0.15 M. Additional analyses included generating LigPlot images before and after simulation, calculating interaction energies, measuring distances between the ligand and catalytic site residues as well as specific residues, computing the root mean square deviation (RMSD) of the protein’s alpha carbon atoms, plotting RMSD distributions, calculating root mean square fluctuations (RMSF) of the alpha carbons, and other relevant assessments33.
RESULT:
Phytochemical constituents of ginger essential oil:
The ginger essential oil contained a total of 40 phytochemicals, as shown Figure 1 and Table 1. These phytochemical components were categorized into two groups of secondary metabolites: terpenoids.
Binding affinity and Interaction of ginger bioactive compounds to orexin receptors:
The selective binding of ginger’s bioactive compounds to orexin receptors was evaluated based on their binding affinity scores. These results were compared with those of lemborexant and suvorexant, which served as control compounds. The findings indicated that ginger bioactive compounds exhibited mild selectivity in binding to orexin receptors (figure 2). Generally, the binding affinities of the ginger compounds were higher than those of the controls. Lemborexant and suvorexant showed selective binding to orexin receptors with low binding affinity values, reflecting strong ligand-receptor interactions. According to the binding affinity scores, suvorexant demonstrated a preference for OX1R (-7.112 kcal/mol), while lemborexant was more selective for OX2R (-10.35 kcal/mol). Among the 40 ginger bioactive compounds tested, alpha-zingiberene exhibited the highest selectivity toward orexin receptors, as indicated by a selectivity score of -0.2342. Its binding affinities were -4.231 kcal/mol for OX1R and -7.255 kcal/mol for OX2R. Notably, alpha-zingiberene showed the strongest binding to OX2R among all ginger compounds, having the lowest binding affinity value for this receptor. Additionally, 6-shogaol, another compound found in ginger essential oil, displayed the strongest interaction with OX1R, with the lowest binding affinity (-5.614 kcal/mol) among all identified ginger bioactive constituents.
Figure 1: TIC Ginger Volatile Oil Compounds using GC-MS QP2010 Plus analysis
Table 1: GC-MS analysis of Ginger Volatile Oil Compounds
|
No. |
Peak number |
RT (min) |
Composition (%) |
Code |
Compound |
Metabolit Secondary |
|
1 |
27 |
24.582 |
18.76 |
CID521253 |
I-zingiberene |
Sesquiterpenoid |
|
2 |
26 |
24.279 |
17.52 |
CID92139 |
alpha-curcumene |
Sesquiterpenoid |
|
3 |
31 |
25.236 |
12.92 |
CID92776 |
alpha-zingiberene |
Sesquiterpenoid |
|
4 |
29 |
24.883 |
7.59 |
CID10104370 |
beta.-bisabolene |
Sesquiterpenoid |
|
5 |
28 |
24.806 |
6.48 |
CID5281516 |
alpha.-Farnesene |
Sesquiterpenoid |
|
6 |
15 |
18.003 |
5.49 |
CID638011 |
citral/ geranial |
Monoterpenoid |
|
7 |
2 |
8.168 |
4.92 |
CID6616 |
camphene |
Monoterpenoid |
|
8 |
13 |
17.143 |
3.92 |
CID643779 |
neral/ cis-citral |
Monoterpenoid |
|
9 |
9 |
14.935 |
2.91 |
CID64685 |
borneol |
Monoterpenoid |
|
10 |
5 |
10.611 |
2.41 |
CID12304570 |
sylvestrene |
Monoterpenoid |
|
11 |
6 |
10.693 |
2.34 |
CID2758 |
eucalyptol/ 1,8-cineol |
Monoterpenoid |
|
12 |
1 |
7.726 |
1.03 |
CID6654 |
alpha-pinene |
Monoterpenoid |
|
13 |
14 |
17.498 |
0.95 |
CID637566 |
geraniol |
Monoterpenoid |
|
14 |
35 |
26.979 |
0.92 |
CID13213649 |
zingiberenol |
Sesquiterpenoid |
|
15 |
19 |
21.197 |
0.78 |
CID12303902 |
copaene |
Sesquiterpenoid |
|
16 |
10 |
15.153 |
0.73 |
CID35725 |
tricyclo[4.3.1.1(3,8)]undecan-1-ol |
Monoterpenoid |
|
17 |
3 |
9.014 |
0.73 |
CID14896 |
beta-pinene |
Monoterpenoid |
|
18 |
40 |
28.351 |
0.67 |
CID5365823 |
dehydronerolidol |
Sesquiterpenoid |
|
19 |
34 |
26.559 |
0.66 |
CID442793 |
6-gingerol |
Sesquiterpenoid |
|
20 |
4 |
9.406 |
0.64 |
CID31253 |
beta.-myrcene |
Monoterpenoid |
|
21 |
20 |
21.674 |
0.61 |
CID6918391 |
beta-elemene |
Sesquitepenoid |
|
22 |
11 |
15.667 |
0.61 |
CID17100 |
alpha-terpineol |
Monoterpenoid |
|
23 |
39 |
28.215 |
0.58 |
CID94378 |
6-paradol |
Monoterpenoid |
|
24 |
8 |
12.793 |
0.57 |
CID6549 |
linalool |
Monoterpenoid |
|
25 |
12 |
16.713 |
0.56 |
CID8842 |
citronellol |
Monoterpenoid |
|
26 |
38 |
27.749 |
0.47 |
CID165258 |
cryptomeridiol |
Sesquitepenoid |
|
27 |
23 |
23.531 |
0.45 |
CID5317319 |
cis-beta.-farnesene |
Sesquitepenoid |
|
28 |
16 |
18.509 |
0.41 |
CID6448 |
bornyl acetate |
Monoterpenoid |
|
29 |
37 |
27.276 |
0.41 |
CID6428444 |
trans-sesquisabinene hydrate |
Sesquitepenoid |
|
30 |
33 |
25.839 |
0.40 |
CID5281794 |
6-shogaol |
Monoterpenoid |
|
31 |
17 |
18.619 |
0.39 |
CID8163 |
2-undecanone |
Monoterpenoid |
|
32 |
22 |
23.056 |
0.35 |
CID15094 |
gamma-muurolene |
Sesquitepenoid |
|
33 |
32 |
25.405 |
0.35 |
CID5352437 |
(E)-gamma-bisabolene |
Sesquitepenoid |
|
34 |
24 |
23.596 |
0.34 |
CID12315492 |
beta-sesquiphellandrene |
Sesquitepenoid |
|
35 |
7 |
12.725 |
0.32 |
CID84825 |
rosefuran/ alpha.-naginatene |
Monoterpenoid |
|
36 |
21 |
22.052 |
0.28 |
CID25147318 |
sesquithujene |
Sesquitepenoid |
|
37 |
30 |
25.059 |
0.28 |
CID91747213 |
cis-muurola-4(15),5-diene |
Sesquitepenoid |
|
38 |
18 |
20.900 |
0.25 |
CID519960 |
(+)-cyclosativene/ cyclosativene |
Sesquitepenoid |
|
39 |
36 |
27.192 |
0.19 |
CID13970960 |
(+/-)-dihydro-ar-turmerone |
Sesquitepenoid |
|
40 |
25 |
23.851 |
0.18 |
CID10899740 |
alloaromadendrene |
Sesquitepenoid |
Figure 2: Selectivities of Binding Affinities Compounds of Ginger Volatile Oil toward Orexin Receptor (OX-1/ OX-2)
Based on binding affinity and selectivity, alpha-zingiberene emerged as the most potent inhibitor of orexin receptors (OXs). The interactions between alpha-zingiberene and the OX1R and OX2R receptors were further examined through 3D visualization (figure 3[c][d]). The binding of alpha-zingiberene to OX1R involved interactions with the amino acid residues Lys1190 and Arg1191 (figure 3[c]) compared to the control compounds suvorexant and lemborexant, alpha-zingiberene exhibited weaker binding, as it engaged fewer amino acid residues in the complex. These results suggest that alpha-zingiberene targets OX1R via a different binding mode than suvorexant and lemborexant (figure 3[a]-[b]), interacting with distinct amino acid residues.
Figure 3: Interaction among suvorexant, lemborexant, and alpha-zingiberene with orexin 1/2 receptor, (a) 4zj8_the human OX-1 receptor with suvorexant (SOV); (b) 7xrr_the human OX2 receptor with lemborexant; (c) 4zj8_the human OX1 receptor with alpha-zingiberene (CID 92776); and (d) 7xrr_the human OX2 receptor with alpha-zingiberene (CID 92776)
Molecular dynamic alpha-zingiberene-OX-2:
Molecular dynamics simulation analysis is essential because it provides a deep understanding of molecular behavior, aids in material prediction and engineering, and complements and explains experimental data that cannot be obtained by conventional methods alone. A demonstration of molecular dynamics simulation using the alpha-zingiberene-OX2R complex is shown in Figure 4.
DISCUSSION:
According to the GC-MS analysis (figure 1 and table1), I-zingiberene and alpha-zingiberene were identified as the most abundant bioactive compounds in the ginger essential oil, comprising 18.76% and 12.92% of the composition, respectively.
Orexins, comprising OX1 and OX2, are neuropeptides that regulate feeding behavior and various neurological functions. Their mechanisms influence the sleep-wake cycle, neuroendocrine activities, glucose metabolism, energy balance, and stress adaptation7,34. Therefore, targeting orexin receptors to modulate the sleep-wake cycle is considered highly promising13,19,35. Suvorexant and emborexant are widely used dual orexin receptor antagonists for treating insomnia8,36. This study, however, demonstrates that the natural bioactive compound alpha-zingiberene from ginger exhibits comparable potency to these commercial orexin receptor antagonists.
[b] RMSD 7xrr_orexin-2-alpha-zingiberene
[c] RMSF 7xrr_orexin-2-alpha-zingeberene
Figure 4: [a] Compare simulation of molecular dynamic 4zj8_orexin-1-alpha-zingiberene and 7xrr_orexin-2-alpha-zingiberene; [b] RMSD parameter of 7xrr/ 4zj8-alpha-zingiberene complexed; and [c] RMSF parameter of 7xrr/ 4zj8-alpha-zingiberene complexed
The identification of potential orexin inhibitors was conducted through in silico analysis by comparing the activity of alpha-zingiberene with that of suvorexant and lemborexant against OX1R and OX2R receptors. Virtual screening via in silico methods serves as a dynamic approach for discovering novel therapeutic agents10,22. The effectiveness of ginger’s bioactive compounds against orexin receptors was evaluated based on their binding affinities and binding sites37.
The interactions of alpha-zingiberene with OX1R and OX2R were further examined through 3D visualization (figure 3[c][d]). The interaction between OX1R and alpha-zingiberene is facilitated by alkyl bonds with Lys A:1190 and Lys A:1191, as well as van der Waals forces involving residues ThrA:1151, AsnA:1151, GluA:1153, ThrA:1154, and AsnA:1187 (figure 3[a]-[d]), compared to suvorexant and lemborexant as controls, alpha-zingiberene forms interaction positions with a distinctly different set of amino acid residues in the complex. These results indicate that alpha-zingiberene targets the OX1 receptor in a manner distinct from suvorexant (figure 3[a]). Alpha-zingiberene also interacts more favorably with OX2R through amino acid residues (figure 3[d]): Pi-sigma (HisA:350), Pi-alkyl (PheA:227), Alkyl (ValA:138, TyrA:317, TyrA:354), and van der Waals (ThrA:111, GlnA:134, SerA:321), which are similar to those involved in the interaction with lemborexant (figure 3[b]: conventional hydrogen bond: GlnA:134; Halogen (Fluorine): Ile320, ArgA:1016; Pi-Pi T-shaped: Phe227; Pi-sigma: IleA:320; and van der Waals: ThrA:111, ValA:114, IleA:130, TrpA:120, ProA:131, ThrA:135, SerA:321, AsnA:324, HisA:350). Similarly, other ginger volatile oils, such as dehydronerolidol, alpha-curcumene, alpha-copaene, alpha-pinene, alpha-terpineol, beta-bisabolene, beta-citronellol, geraniol, linalool, and others, also tend to inhibit the human orexin-2 receptor, although their binding affinities are lower than that of alpha-zingiberene, as shown in Figure 2. Therefore, alpha-zingiberene is predicted to be the best ginger oil compound for inhibiting the human orexin-2 receptor.This study found that alpha-zingiberene exhibited the highest selective binding activity toward both OX1R and OX2R among all identified ginger compounds. As an orexin antagonist, alpha-zingiberene likely targets OX1R and OX2R with binding behaviors moderately similar to those of lemborexant. alpha-zingiberene may inhibit the interaction of orexins with their respective G-protein-coupled receptors (GPCRs), specifically OX1R with Hcrtr-1 and OX2R with Hcrtr-2. An orexin antagonist, such as a Dual Orexin Receptor Antagonist (DORA), binds to these receptors without activating them, effectively blocking orexin from binding and initiating its downstream signaling cascade. By acting as an orexin antagonist, it can block the activation of phospholipase C (PLC), phospholipase A (PLA), phospholipase D (PLD), and adenylyl cyclase (AC), leading to reduced cytosolic Ca2+ levels and suppression of downstream signaling cascades7.
Molecular dynamics simulation analysis is essential because it provides a deep understanding of molecular behavior, aids in material prediction and engineering, and complements and explains experimental data that cannot be obtained by conventional methods alone. A demonstration of molecular dynamics simulation using the alpha-zingiberene-OX2R complex is shown in Figure 4. The stability of the protein–ligand complexes was assessed by calculating the Root Mean Square Deviation (RMSD) and the Root Mean Square Fluctuation (RMSF) of the protein backbone from its initial to final conformation. RMSD (Root Mean Square Deviation) is a fundamental metric in molecular dynamics (MD) simulations used to quantify the average structural deviation of a molecular system, such as a protein or nucleic acid over time relative to a reference structure, often the initial or an experimentally determined conformation. RMSD simulations showed that the alpha-zingiberene-OX2R complex reached overall stability after 2 ns (nanoseconds), with an average stable RMSD of 1.9 Å (figure 4[b]). A lower RMSD value at 2 ns means that the molecule (protein, ligand, or their complex) maintains its initial conformation well and is more dynamically stable. A stable and low RMSD value (around 0.1–0.2 nm or 1–2 Å) during the simulation is a key indicator of molecular stability in molecular dynamics simulations. These results indicate that the alpha-zingiberene-OX2R complex is more stable during the simulation time than other ginger oil compounds. The Root Mean Square Fluctuation (RMSF) measures the flexibility of protein residues by quantifying their deviation from average positions during the simulation. The root mean square fluctuation (RMSF) value is the fluctuation of protein residues from their mean positions during the simulation, which describes the flexibility of the protein structure. Based on the timeline results, alpha-zingiberene interacts effectively with the orexin-2 receptor binding site (figure [4c]. alpha-zingiberene interacts with the active site residues and the A-domain site, and reduces the site fluctuation in orexin-2R. Furthermore, Figure 4a illustrates the timeline interaction of alpha-zingiberene with the active site of orexin-2 at simulation time.
The binding energy of alpha-zingiberene-orexin-2R complex using MMPBSA and MMGBSA solvent model was reported by Open MM on Google Colab. The MMPBSA binding energy (∆G) of alpha-Zingiberene-OX2R complex is 0.0226 ± 3.1674 kcal/ mol and the MMGBSA binding energy (∆G) of alpha-Zingiberene-OX2R complex is -8.9200 ± 2.3693 kcal/ mol. The difference in the binding energy value between MMPBSA and MMGBSA could be due to polar and nonpolar contributions to ΔG, ΔT and ΔS, and type of bonds, angles, dihedral energies, electrostatic energies and van der Waals energies. MMGBSA/PBSA binding energy (∆G) is an estimate of the free energy of binding between ligand (alpha-zingiberene) and receptor (OX1R or OX2R) calculated using the MMGBSA/MMPBS method. A more negative ∆G value indicates a stronger and more stable bond between the ligand and the receptor.
The ∆G value of OX2R is negative, indicating that alpha-zingiberene forms a stronger and more stable bond to OX2R. Several key contributing factors are the characteristics of the binding pocket, the geometric and chemical compatibility between the ligand and the receptor, and the dynamics and flexibility of the protein. A look at the structure of the binding pocket, illustrated in Figure 4[a], shows that critical residues in OX1R (e.g., Val160, Thr161, Arg162, Phe163, Arg164, Ala127 initially, and Leu4, Arg5, Trp8, His13 during the simulation) play a crucial role in stabilizing interactions with ligands such as alpha-zingiberene. These residues engage in hydrophobic interactions and sometimes form hydrogen bonds favorable for certain ligands. Analysis of the geometric and chemical compatibility between the ligand and the receptor indicates that alpha-zingiberene fits better spatially and chemically within the binding pocket of OX2R, resulting in lower (more negative) binding energy and thus greater stability. Conversely, a geometric or polarity mismatch in the OX2R pocket will weaken the interaction, leading to a less negative ∆G value. Furthermore, the dynamics and flexibility of the OX2R protein affect the stability of the ligand-receptor complex during the simulation, as reflected in the MMGBSA binding energy values and their deviations. The standard deviation (±) values indicate that OX2R exhibits a relatively low deviation indicating a more consistent but generally weaker interaction. Overall, the binding energy of the alpha-zingiberene–OX2R complex indicates stability and is consistent with the RMSD, RMSF, and binding affinity analyses38.
Ginger is a famous medicinal herb known for its biological activities. Bioactive compounds of ginger were proven for the activities in neurophysiological mechanisms and in the blood-brain barrier18,19,39. Gingerols and shogaols, as the major constituents in ginger, performed activities including migraine relief, neuroprotective, monoamine oxidase A (MAO-A), and GABA modulator40. Previously, treatment of Z. officinale hydroethanolic extract successfully exerted antiseizure activity in the PTZ-kindling mice model41. Treatment with ginger also effectively upregulated GABA level in the hippocampus and cortex of the senile female rats19,42. This study identified the potency of alpha-zingiberene to act as orexin antagonist via modulation of GABA, similarly as lemborexant36. Treatment with alpha-zingiberene via OX1R/OX2R may increase the sleep drive and inhibit wake-promoting neurons, improving total amount of sleep and daytime alertness35. Ginger, as natural herb, may exhibit less adverse effects than drugs13,11 43. However, the physiological effect of alpha-zingiberene remains unknown, thus it required further studies to evaluate the potency.
CONCLUSION:
This study identified novel orexin antagonists among the natural bioactive compounds from ginger essential oil. Molecular docking results indicated that these compounds were generally selective for the OX2R receptor. Molecular dynamics simulations also showed that alpha-zingiberene formed a stronger and more stable interaction with OX2R. These findings demonstrate the potential of ginger essential oil compounds, particularly alpha-zingiberene, as promising therapeutic agents for the treatment of insomnia.
CONFLICT OF INTEREST:
The authors have no conflicts of interest regarding this investigation.
ACKNOWLEDGMENTS:
The authors acknowledge to Ministry of Education and Culture for research funding through a national research grant of Indonesia year 2024.
REFERENCES:
1. Burnett A, DiTunnariello N, DeGreeff BL. “I’m on a Rollercoaster”: Women’s Social Construction of Time. Communication Studies. 2020; 71: 148–66. https://doi.org/ 10.1080/10510974.2019.1702073.
2. Lippke S, Schalk TM, Kühnen U, Shang B. Pace of life and perceived stress in international students. PsyCh Journal. 2021; 10: 425–36. https://doi.org/10.1002/pchj.426.
3. Maanickha Chelvi KS, Karthik Nagarajan S, Kanagarajan A. Insomnia Management Clinical Trials by various Conventional Therapies – A Review. Research Journal of Pharmacy and Technology. 2023; 16(12): 6071-6076. https://doi.org/ 10.52711/0974-360X.2023.00985
4. Suneesh PM, Thangavel K, Doss KJJ. A Study to Assess the Effectiveness of Warm Foot Bath on Sleep Onset Time among Cancer Patients with Insomnia in Nath Lal Parekh Cancer Hospital at Rajkot. Asian J Nurs Educ Res. 2017; 7(2): 231-234. https://doi.org/10.5958/2349-2996.2017.00048.9.
5. Sharma M, Prakash K, Priya J. Effectiveness of “Warm Footbath” on fatigue and insomnia in patients undergoing radiotherapy. Int J Nurs Educ Res. 2016; 3(4): 363-366. https://10.5958/2454-2660.2016.00064.8
6. Chieffi S, Carotenuto M, Monda V, Valenzano A, Villano I, Precenzano F. Orexin System: The Key for a Healthy Life. Front Physiol. 2017; 8: 357. https://doi.org/10.3389/fphys.2017.00357.
7. Wang C, Wang Q, Ji B, Pan Y, Xu C, Cheng B, et al. The Orexin/Receptor System: Molecular Mechanism and Therapeutic Potential for Neurological Diseases. Front Mol Neurosci. 2018; 11: 220. https://doi.org/10.3389/fnmol.2018.00220.
8. Khazaie H, Sadeghi M, Khazaie S, Hirshkowitz M, Sharafkhaneh A. Dual orexin receptor antagonists for treatment of insomnia: A systematic review and meta-analysis on randomized, double-blind, placebo-controlled trials of suvorexant and lemborexant. Front Psychiatry. 2022; 13: 1070522. https://doi.org/10.3389/fpsyt.2022.1070522.
9. Muehlan C, Roch C, Vaillant C, Dingemanse J. The orexin story and orexin receptor antagonists for the treatment of insomnia. Journal of Sleep Research. 2023; 32: e13902. https://doi.org/10.1111/jsr.13902.
10. Rashmi KS, Maynnavar SK, Arunkumar N, Bhagyalakshmi K, Poojary D, Kadur K, et al. The regulatory role of neurotransmitters in orexin B induced nucleus accumbens dependent alcohol ingestion in wistar albino rats. Res J Pharm Technol. 2020; (13)12: 6224-6230. https://doi.org/ 10.5958/0974-360X.2020.01085.9.
11. Ardianto C, Judoko AL, Aprilia M, Ratri DMN, Ariyani T, Rahmadi M, et al. Effect of Andrographolide and Resveratrol on OX1R and Prepro-orexin mRNA expression in CIPN-induced hypothalamus of mice with oxaliplatin. Res J Pharm Technol. 2022; 15(10): 4765-4771. https://doi.org/ 10.52711/0974-360X.2022.00800
12. Mogavero MP, Silvani A, Lanza G, DelRosso LM, Ferini-Strambi L, Ferri R. Targeting Orexin Receptors for the Treatment of Insomnia: From Physiological Mechanisms to Current Clinical Evidence and Recommendations. NSS. 2023; 15: 17–38. https://doi.org/10.2147/NSS.S201994.
13. Jasim RF, Allwsh TA. Orexin a hormone and its relation to coronary heart diseases. Res J Pharm Technol. 2021; 14(3): 1417-1422. https://doi.org/10.5958/0974-360X.2021.00253.5.
14. Ajayi BO, Olajide TA, Olayinka ET. 6-gingerol attenuates pulmonary inflammation and oxidative stress in mice model of house dust mite-induced asthma. Advances in Redox Research. 2022; 5: 100036. https://doi.org/10.1016/j.arres.2022.100036.
15. Arcusa R, Villaño D, Marhuenda J, Cano M, Cerdà B, Zafrilla P. Potential Role of Ginger (Zingiber officinale Roscoe) in the Prevention of Neurodegenerative Diseases. Front Nutr. 2022; 9: 809621. https://doi.org/10.3389/fnut.2022.809621.
16. Bare Y, Maulidi A, Sari DRT, Tiring SSND. Studi in Silico Prediksi Potensi 6-Gingerol sebagai inhibitor c-Jun N-terminal kinases (JNK). JJMS. 2019; 1: 59–63. https://doi.org/10.36873/jjms.v1i2.211.
17. Gawel K, Kukula-Koch W, Banono NS, Nieoczym D, Targowska-Duda KM, Czernicka L, et al. 6-Gingerol, a Major Constituent of Zingiber officinale Rhizoma, Exerts Anticonvulsant Activity in the Pentylenetetrazole-Induced Seizure Model in Larval Zebrafish. IJMS. 2021; 22: 7745. https://doi.org/10.3390/ijms22147745.
18. Simon A, Darcsi A, Kéry Á, Riethmüller E. Blood-brain barrier permeability study of ginger constituents. Journal of Pharmaceutical and Biomedical Analysis. 2020; 177: 112820. https://doi.org/10.1016/j.jpba.2019.112820.
19. Wilson SA, Nagaraj A, Vaidyanathan L. Calocybe indicia extract modifies gaba and serotonin levels to alleviate anxiety in experimental adult zebrafish models. Res J Pharm Technol. 2021; 14(9): 4789-4794. https://doi.org/ 10.52711/0974-360X.2021.00833
20. Sawhney SK, Narayan C, Mishra A, Singh M, Kaur A. Molecular Docking, ADME and Toxicity Study of Dibenzo-α-pyrone derivatives for GABA and NMDA receptors for their antiepileptic activity. Res J Pharm Technol. 2024; 17(1): 340-346. https://doi.org/ 10.52711/0974-360X.2024.00053
21. Uddin J, Ahmed H, Asiri YI, Kamal GM, Musharraf SG. Ginger essential oil: Chemical composition, extraction, characterization, pharmacological activities, and applications. Essential Oils, Elsevier. 2023; 345–76. https://doi.org/10.1016/B978-0-323-91740-7.00014-1.
22. Hellmann J, Drabek M, Yin J, Gunera J, Pröll T, Kraus F, et al. Structure-based development of a subtype-selective orexin 1 receptor antagonist. Proc Natl Acad Sci USA. 2020; 117(30): 18059-18067. https//doi.org/ 10.1073/pnas.2002704117.
23. Asada H, Im D, Hotta Y, Yasuda S, Murata T, Suno R, et al. Molecular basis for anti-insomnia drug design from structure of lemborexant-bound orexin 2 receptor. Structure. 2022; 30: 1582-1589.e4. https://doi.org/10.1016/j.str.2022.11.001.
24. Trott O, Olson AJ. AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J Comput Chem. 2010; 31: 455–61. https://doi.org/10.1002/jcc.21334.
25. Liao A, Li L, Wang T, Lu A, Wang Z, Wang Q. Erratum: Discovery of Phytoalexin Camalexin and Its Derivatives as Novel Antiviral and Antiphytopathogenic-Fungus Agents. J. Agric. Food Chem. 2022; 70: 2554-2563. https://doi.org/10.1021/acs.jafc.1c07805).
26. Sari D, Safitri A, Cairns J, Fatchiyah F. Virtual screening of black rice anthocyanins as antiobesity through inhibiting TLR4 and JNK pathway. J Phys: Conf Ser. 2020; 1665: 012024. https://doi.org/10.1088/1742-6596/1665/1/012024.
27. Raharjo SJ, Kikuchi T. Molecular dynamic screening sesquiterpenoid pogostemon herba as suggested cyclooxygenase inhibitor. Acta Inform Medica. 2016; 24(5): 332–7. https://doi.org/ 10.5455/aim.2016.24.332-337.
28. Ananda AN, Triawanti T, Setiawan B, Makati AC, Putri JA, Raharjo SJ. In silico study of the flavonoid compound of Sauropus androgynus leaves ON RNA-Dependent RNA polymerase (RdRp) SARS-CoV-2. Asp Mol Med. 2024; 3(6); 1-8. https://doi.org/10.1016/j.amolm.2023.100032
29. Utami PD, Setianingsih H, Sari DRT. In Silico Approach Triterpene Glycoside of H. atra Targeting Orotidine 5-Monophosphate Decarboxylase Protein (PfOMPDC) in P. falciparum Infection Mechanism. BioMed Research International 2024; 2024: 1–14. https://doi.org/10.1155/2024/5924799.
30. Dodda LS, Vilseck JZ, Tirado-Rives J, Jorgensen WL. 1.14*CM1A-LBCC: Localized Bond-Charge Corrected CM1A Charges for Condensed-Phase Simulations. J Phys Chem B. 2017; 12(15): 3864–70.https//doi.org/10.1021/acs.jpcb.7b00272.
31. Dodda LS, Cabeza de Vaca I, Tirado-Rives J, Jorgensen WL. LigParGen web server: an automatic OPLS-AA parameter generator for organic ligands. Nucleic Acids Research. 2017; 45: W331–6. https://doi.org/10.1093/nar/gkx312.
32. Jorgensen WL, Tirado-Rives J. Potential energy functions for atomic-level simulations of water and organic and biomolecular systems. Proc Natl Acad Sci. 2005; 102: 6665–70.
33. Arantes PR, Polêto MD, Pedebos C, Ligabue-Braun R. Making it Rain: Cloud-Based Molecular Simulations for Everyone. J Chem Inf Model. 2021; 61: 4852–6. https://doi.org/10.1021/acs.jcim.1c00998.
34. Cao M, Chow M. The hypocretin/orexin system in sleep disorders: preclinical insights and clinical progress. NSS. 2016: 81. https://doi.org/10.2147/NSS.S76711.
35. Yamanaka M, Iwata H, Masuda K, Araki M, Okuno Y, Okamura M, et al. A novel orexin antagonist from a natural plant was discovered using zebrafish behavioural analysis. European Review for Medical and Pharmacological Sciences 2020; 24: 5127–39. https://doi.org/ 10.26355/eurrev_202005_21207.
36. Reddy PPK, Hiremath JS, Nyamagoud SB, Swamy AHV. A comprehensive review of the use of lemborexant in treating insomnia. MGMJ. 2024; 11: 146–51. https://doi.org/10.4103/mgmj.mgmj_286_23.
37. Krisnamurti GC, Fatchiyah F. The Biological Function Prediction of The 10-gingerol Compound of Ginger in Inhibiting Cyclooxygenase-2 Activity. J Pure App Chem. Res. 2020; 9: 222–32. https://doi.org/10.21776/ub.jpacr.2020.009.03.547.
38. Genheden S and Ryde U. The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities. Expert Opin Drug Discov. 2015; 10(5): 449–461. doi: 10.1517/17460441.2015.1032936
39. Sharma R, Kumari N, Ashawat M., Verma CP. Standardization and phytochemical screening analysis for herbal extracts: Zingiber officinalis, Rosc., Curcuma longa Linn., Cinnamonum zeylanicum Nees., Piper longum, Linn., Boerhaavia diffussa Linn. Asian J Pharm Technol. 2020; 14(1): 341-346. https://doi.org/ 10.52711/0974-360X.2024.00053
40. Kukula-Koch W, Koch W, Czernicka L, Głowniak K, Asakawa Y, Umeyama A, et al. MAO-A Inhibitory Potential of Terpene Constituents from Ginger Rhizomes. A Bioactivity Guided Fractionation. Molecules. 2018; 23: 1301. https://doi.org/10.3390/molecules23061301.
41. Hosseini A, Mirazi N. Alteration of pentylenetetrazole-induced seizure threshold by chronic administration of ginger (Zingiber officinale) extract in male mice. Pharmaceutical Biology. 2015; 53: 752–7. https://doi.org/10.3109/13880209.2014.942789.
42. Shukla SS, Saraf S, Saraf S. Development and validation of spectrophotometric fingerprint method of 6-gingerol in herbal formulation: Talisadi churna. Res J Pharm Technol. 2012; 5(1): 1391-40.
43. Asma ST, Acaroz U, Imre K, Morar A, Shah SRA, Hussain SZ, et al. Natural Products/ Bioactive Compounds as a Source of Anticancer Drugs. Cancers. 2022; 14: 6203. https://doi.org/10.3390/cancers14246203.
|
Received on 13.07.2025 Revised on 24.11.2025 Accepted on 11.02.2026 Published on 01.07.2026 Available online from July 04, 2026 Research J. Pharmacy and Technology. 2026;19(7):3325-3332. DOI: 10.52711/0974-360X.2026.00473 © RJPT All right reserved
|
|
|
This work is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License. Creative Commons License. |
|